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1. HISTORICAL INTRODUCTION 



While the Earth's magnetism has been studied for 
centuries starting with the first scientific monograph of 
Gilbert (1600), the question of the magnetism of other 
planets had received scant attention until recently be- 
cause of the lack of relevant observations. Only in 1955 
clear evidence for the existence of a planetary mag- 
netic field other than the geomagnetic one was obtained 
through the observation of the Jovian decametric radio 
waves (Burke and Franklin, 1955). Since it had been 
more or less accepted until the end of the 19-th century 
that geomagnetism arises from the remnant magnetiza- 
tion of the Earth similar properties may have been as- 
sumed for the other terrestrial planets and the Moon. 
This view lost its appeal, however, when it became evi- 
dent that the Curie-temperature is exceeded in the Earth 
below a depth of about 30 km. Ferromagnetic materials 
in the Earth's crust could thus explain only magnetic 
fields with a relatively short wavelength. 

The current period of intense research on the mag- 
netism of planets other than that of the Earth started 
with the first detailed measurement of Jupiter's mag- 
netic field by the Pioneer 10 space probe in 1973 and 
the discovery of Mercury's magnetism by Mariner 10 in 
1974. In the early seventies also the development of the 
theory of magnetohydrodynamic dynamos had started in 
which the reaction of the Lorentz force of the generated 
magnetic field is taken into account in physically real- 
istic configurations (Childress and Soward, 1972; Busse, 
1973; Soward, 1974). Until that time dynamo theoreti- 
cians had focused their attention on the kinematic prob- 
lem in which the possibility of growing magnetic fields 
driven by somewhat arbitrarily chosen velocity fields is 
considered. It must be remembered that only a few years 
earlier it had been demonstrated by Backus (1958) and 



Herzenberg (1958) in a mathematically convincing way 
that the homogeneous dynamo process of the generation 
of magnetic fields in a singly connected domain of an 
electrically conducting fiuid is indeed possible. Doubts 
about the feasibility of this process which had first been 
proposed by Larmor (1919) had persisted after Cowling 
(1934) had proved that purely axisymmetric fields could 
not be generated in this way. 

The complexity of the magnetohydrodynamic dynamo 
problem described by the nonlinearly coupled Navier- 
Stokes equations and the equation of magnetic induction 
had prevented progress in understanding planetary dy- 
namos through analytical solutions. Only the advent of 
powerful enough computers in the 1990-ies has allowed 
to solve numerically the coupled three-dimensional par- 
tial differential equations through forward integration in 
time. Even today and for the foreseeable future the lim- 
its of computer capacity will permit the exploration of 
only a small fraction of the parameter space of interest 
for the understanding of planetary dynamos. 

In view of the difficulties of a rigorous theory of plane- 
tary dynamos, many attempts have been made to obtain 
simple similarity relationships which would fit the ob- 
served planetary magnetic moments as function of cer- 
tain properties of the planets. Some early proponents 
have gone as far as claiming the existence of a "magnetic 
Bode's law" corresponding to a relationship between the 
magnetic moment and size or angular momentum of a 
planet in analogy to the Titius-Bode law for the radii 
of the orbits of the planets. Just as in the latter case, 
however, attempts to derive a magnetic Bode's law from 
basic physical principles have failed. 

Other proposals have taken into account physical 
forces. Since a common ingredient of planetary dynamos 
is the existence of a fiuid part of the core with a suffi- 
ciently high electrical conductivity the latter parameter 
together with the core radius and the angular velocity of 
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the planetary rotation usually enter the similarity rela- 
tionships such as those proposed by Hide (1974), Busse 
(1976), Jacobs (1979) and Dolginov (1977). Malkus 
(1968, 1996) has argued for the precession of the Earth 
as the cause of geomagnetism and he and Vanyo (1984) 
have demonstrated through laboratory experiments that 
precession and tides may cause turbulent motions in fluid 
planetary cores. Dolginov (1977) proposed a scaling law 
for the precessional origin of all planetary magnetic fields. 
While a dynamo driven by turbulent flows caused by pre- 
cession and tides can not be easily excluded in the case 
of the Earth (Tilgner, 2005), it is much less likely in the 
case of other planets such as Uranus for which preces- 
sional torques are rather minute. Just as a common pre- 
cessional origin of planetary magnetism is not regarded 
as feasible, so have all other proposed similarity relation- 
ships lost in appeal and are no longer seriously consid- 
ered. We shall return, however, to scaling relationships 
based more directly on the basic equations in section 6. 



2. GENERAL REMARKS ON THE DYNAMO 
THEORY OF PLANETARY MAGNETISM 

Since the proposal of the gcodynamo as the cause of the 
Earth's magnetism had been in doubt for a long time be- 
fore 1958, numerous alternative proposals had been made 
in the literature. Among these only the possibility that 
thermoelectric currents may generate a planetary mag- 
netic field is still discussed in the case of Mercury (Steven- 
son, 1987; Giampieri and Balogh, 2002). For a discussion 
of the failings of the various proposals for non-dynamo 
origins of planetary magnetic fields we refer to the papers 
of BuUard (1949) and Stevenson (1983). Although the 
dynamo hypothesis of the origin of planetary magnetism 
is not without difiiculties, it is the only one considered 
seriously at the present time with the possible exception 
of the just mentioned case of Mercury. 

Dynamos generally convert mechanical energy into 
magnetic one. In contrast to the technical dynamo which 
is characterized by a multiply connected region of high 
electrical conductivity, i.e. it depends on an appropriate 
wiring, planetary dynamos arc referred to as homoge- 
neous dynamos since they operate in a singly connected 
domain of high electrical conductivity. Since fiows in 
planetary cores with active dynamos arc usually turbu- 
lent the small scale structure of the magnetic field is cor- 
respondingly chaotic. The large scale structure, however, 
can be quite regular. One distinguishes "steady" and os- 
cillatory dynamos. The most famous example of the lat- 
ter kind is the solar dynamo which exhibits a well defined 
period of about 22 years. The gcodynamo, on the other 
hand, is a "steady" dynamo, even though it varies in its 
amplitude by a factor of two or more on the magnetic 
diffusion time scale and reverses its polarity on a longer 
time scale. A measure of the magnetic diffusion time is 



given by the decay time, td = (r^ir^/ir'^, of the magnetic 
field in the absence of fluid motions. Here a and refer 
to the electrical conductivity and the magnetic perme- 
ability of the planetary core of radius tq. In the case of 
the Earth the decay time is of the order of 20 000 years, 
but it may vary between a few hundred and a million 
years for other examples of planetary dynamos. 

The theory of homogeneous dynamos is based on 
Maxwell's equations for the magnetic flux density B, the 
electric current density j and the electric field E in the 
magnetohydrodynamic approximation in which the dis- 
placement current is neglected. This approximation is 
highly accurate as long as the fluid velocity is small com- 
pared to the velocity of light which is certainly the case 
for all planetary applications. These equations together 
with Ohm's law for a moving conductor are given by 



V-B = 0, 
V X B = H3, 



= -V xE, (la) 
j = (j{vxB + E), (lb) 



where fi is the magnetic permeability of the fluid and a 
is its electrical conductivity. These "pre-Maxwcll" equa- 
tions have the property that they are invariant with re- 
spect to a Galilei transformation, i. e. the equations 
remain unchanged in a new frame of reference moving 
with the constant velocity vector V relative to the orig- 
inal frame of reference. Indicating the variables of the 
new frame by a prime we find 



B' = B, 



dt' 
E' = E 



dt 
V X B, 



3 =3 



(2a) 
(2b) 



This invariance is the basis for the combination in MHD 
of equations (1) with the equations of hydrodynamics 
in their usual non-relativistic form. It is remarkable 
that this invariance does not only hold with respect to 
a Galilei transformation, but with respect to a transfor- 
mation to a rotating frame of reference as well. In that 
case V is replaced by r2 x r in equations (2), but when 
^ is operating on any vector a the term — n x o must 
be added on the right hand side, since even a constant 
vector field becomes time dependent when seen from a 
rotating frame unless it is parallel to ft. 

Elimination of E and j from equation (1) yields the 
equation of magnetic induction 

V X (—V X b] = -^B + V X (v X B) (3) 
\c7fj, J at 

which for a solenoidal velocity field v and a constant mag- 
netic diffusivity A = (cr/i)~^ can be further simplified, 



d_ 

dt 



v-V\B- XV^B ^BVv. 



(4) 



This equation has the form of a heat equation with the 
magnetic field line stretching term on the right hand side 
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FIG. 1: Geometrical configuration of the problem. Apart 

of the outer spficrical surface is removed to expose tlie 
interior of the shell to which the conducting fluid is confined. 



acting as a heat source. This interpretation is especiaUy 
useful for the dynamo problem. In order that a magnetic 
field B may grow, the term on the right hand side of (4) 
must overcome the effect of the magnetic diffusion term 
on the left hand side. Using a typical velocity U and a 
typical length scale d, the ratio of the two terms can be 
estimated by the magnetic Reynolds number Rm, 



Ud/X = Rm. 



(5) 



Only when Rm is of the order one or larger may growing 
magnetic fields become possible. 

In the following we shall first consider the mathemat- 
ical formulation of the problem of convection driven dy- 
namos in rotating spherical shells in a simple form in 
which only the physically most relevant parameters are 
taken into account. Before discussing dynamo solutions 
in section 5 we shall briefly outline in section 4 properties 
of convection it the absence of a magnetic field. Appli- 
cations to various planets and moons will be considered 
in section 6 of this article and some concluding remarks 
are given in section 7. 



outer radius of the spherical shell and d is its thickness. 
AT is the temperature difference between the boundaries 
in the special case /3 = 0. In the case AT = the static 
temperature distribution T5 corresponds to that of a ho- 
mogeneously heated sphere with the heat source density 
proportional to the parameter /3. The gravity field is 
given by g = —jdr where r is the position vector with 
respect to the center of the sphere and r is its length 
measured in units of d. 

In addition to d, the time d'^/v, the temperature 
v^l^Oid^ and the magnetic flux density v(jiq)^^^ j d are 
used as scales for the dimensionless description of the 
problem where v denotes the kinematic viscosity of the 
fluid, K its thermal diffusivity and q its density. The 
Boussinesq approximation is used in that q is assumed to 
be constant except in the gravity term where its temper- 
ature dependence given by a = —{dQ/dT)/Q = const, is 
taken into account. The dimensionless equations of mo- 
tion, the heat equation for the deviation of the temper- 
ature fleld from the static distribution and the equation 
of magnetic induction thus assume the form 

V^v + B ■ VB + re - Vtt 

= P~^{dtV + v-Vv)+T\<ixv (6a) 
V • w = (6b) 

V^e + [Ri + ReW'^O^ - rf)-^] r ■ v 

= P{dt + V ■ V)e (6c) 

Pm + • B - V'B = PmB ■ Vv (6d) 

where k is the unit vector in the direction of the axis 
of rotation and where Vtt includes all terms that can be 
written as gradients. The Rayleigh numbers R4 and Re, 
the Coriolis parameter t, the Prandtl number P and the 
magnetic Prandtl number Pm are defined by 



Ri 



2nd^ 



, Re 



ajATd'^ 



P = - P 

V K 



V 

A 



(7a) 
(7b) 



3. MATHEMATICAL FORMULATION OF THE 
PROBLEM OF SPHERICAL DYNAMOS 

A sketch of the geometrical configuration that will be 

considered is shown in figure 1. For the equations de- 
scribing convection driven dynamos in the frame rotat- 
ing with the angular velocity f2 we use a standard for- 
mulation which has also been used for a dynamo bench- 
mark (Christensen et al. 2001). But we assume different 
scales and assume that a more general static state exists 
with the temperature distribution Tg = Tq — (id^r^ /2 + 
l\Trir~^(l — rf)~^ where rj denotes the ratio of inner to 



For simplicity Re = will be assumed unless indicated 
otherwise. The notation R = Ri will thus be used in the 
following. The Prandtl number P has been added as an 
important parameter of the problem even though P = 
1 is often assumed with the argument that all effective 
diffusivities are equal in turbulent media. The effective 
diffusivities of scalar and vector quantities in turbulent 
fiuid fiow differ in general, however, and largo differences 
in the corresponding molecular diffusivities will not be 
erased entirely in the turbulent case. 

Since the velocity field v as well as the magnetic flux 
density B are solenoidal vector fields, the general repre- 
sentation in terms of poloidal and toroidal components 
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can be used, 

v = V X (V$ X r) + V* X r , (8a) 
B = V X (V/i X r) + Vfif X r . (8b) 

By multiplying the (curl)^ and the curl of equation (6a) 
by r we obtain two equations for $ and ^ 

[(V^ - dt)L2 + Td4,]VH + tQ^ - L^Q 

= -r • V X [V X (v -Wv-B- VB)] (9a) 

[{V^ -dt)L2+Td^]^ -tQ'^ 

= r ■ W X {v ■ Wv - B ■ WB) (9b) 

where dt and denote the partial derivatives with re- 
spect to time t and with respect to the angle of a spher- 
ical system of coordinates r, 6, (j) and where the operators 
1/2 and Q are defined by 

La = -r^V^ + drir'^dr) 

Q = rcoseV^ - {L2 + rdr)(cosedr - r'^ sin Ode) 

The equations for h and g arc obtained through the mul- 
tiplication of equation (6d) and of its curl by r, 

V^La/i = Prn[dtL2h -r-Vx{vxB)] (10a) 
V^L2g = Pm[dtL2g - r • V X (V X (t; X B))] (10b) 

Either rigid boundaries with fixed temperatures as in the 
benchmark case (Christensen et al. 2001), 

$ = dr{r^) = * = e = at r = ri = 77/(1 - r?) 

and at r = ro = (1 — 7?)""^, (11) 
or stress-free boundaries with fixed temperatures, 
$ = = dri'b/r) = 6 = at r = ri = 77/(1 - 77) 

and at r = To = (1 - 77)"^ (12) 

are often used. The latter boundary conditions are as- 
sumed in the following since they allow to cover numer- 
ically a larger region of the parameter space. The case 
77 = 0.4 will be considered unless indicated otherwise, ft 
provides a good compromise for the study of both, the 
regions inside and outside the tangent cylinder. The lat- 
ter is the cylindrical surface touching the inner spher- 
ical boundary at its equator. For the magnetic field 
electrically insulating boundaries are used such that the 
poloidal function h must be matched to the function h^'^'^ 
which describes the potential fields outside the fluid shell 

g = h- h^"'^ = dr{h - h^"'^) = at r = n and r = r,,. 

(13) 



E 


Average density of kinetic energy 


(18) 


Ep, Et 


Energy densities of poloidal 






and toroidal components of motion 


(18) 




Radius ratio of fluid shell 




A 


Elsasser number 


(20) 


M 


Average magnetic energy density 


(19) 


Mp, Mt 


Energy densities of poloidal and toroidal 






components of the magnetic field 


(19) 


Nu 


Nusselt number 


(16) 


P 


Prandtl number 


(7b) 


p 


Magnetic Prandtl number 


(7b) 


R 


Rayleigh number 


(7a) 


Rc 


Critical value of R for onset of convection 




Rm 


Magnetic Reynolds number = PmV2E 




T 


Coriolis parameter 


(7a) 



TABLE I: Important dimensionless parameters. 



But computations for the case of an inner boundary with 
no-slip conditions and an electrical conductivity equal to 
that of the fluid have also been done. The numerical 
integration of equations (2) together with boundary con- 
ditions (4) proceeds with the pseudo-spectral method as 
described by Tilgner and Busse (1997) and Tilgner (1999) 
which is based on an expansion of all dependent variables 
in spherical harmonics for the 9, ^-dependences, i.e. 

$ = ^ (r, t)Pr (cos 6*) exp{im<p} (14) 

and analogous expressions for the other variables, 4", 6, /i 
and g. P;"* denotes the associated Legendre functions. 
For the r-dependence expansions in Chebychev polyno- 
mials are used. 

For the computations to be reported in the following a 
minimum of 33 collocation points in the radial direction 
and spherical harmonics up to the order 64 have been 
used. But in many cases the resolution was increased to 
49 collocation points and spherical harmonics up to the 
order 96 or 128. 

It should be emphasized that the static state v = B = 
Q = represents a solution of equations (6) for all values 
of the Rayleigh numbers Ri and Re, but this solution is 
unstable except for low or negative values of the latter 
parameters. Similarly, there exist solutions with B — 0, 
but f 7^ 0, 7^ 0, for sufficiently large values of either Ri 
or Re or both, but, again, these solutions are imstable 
for sufficiently large values of Pm with respect to distur- 
bances with B 0. Dynamo solutions as all solutions 
with B ^ are called are thus removed by at least two 
bifurcations from the basic static solution of the problem. 

Finally we present in table 1 a list of the most impor- 
tant parameters used in the following sections. 
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4. CONVECTION IN ROTATING SPHERICAL 
SHELLS 

For an introduction to the problem of convection in 
spherical shells we refer to the review of Busse (2002a) . 
Convection tends to set in first outside the tangent cylin- 
der in the form of thermal Rossby waves for which the 
Coriolis force is balanced almost entirely by the pressure 
gradient. The model of the rotating cylindrical annulus 
has been especially useful for the understanding of this 
type of convection. A rough idea of the dependence of 
the critical Rayleigh number Rc for the onset of convec- 
tion on the parameters of the problem can be gained from 
the expressions derived from the annulus model (Busse, 
1970) 

Rc = 3(^-^^y {t&n9m)ir~h-i (15a) 

/ Pt \^ 2 _ 1 

"^'^^iiTpj " ^^^^^ 

where refers to the mean radius of the fluid shell, = 
{ri+ro)/^, and 6m to the corresponding colatitude, 9m = 
arcsin {vmi^ ~ v))- The azimuthal wavenumber of the 
preferred mode is denoted by rric and the corresponding 
angular velocity of the drift of the convection columns in 
the prograde direction is given by Uc/rric. 
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FIG. 2: Critical Rayleigh number Rc (thin solid linos) and 
frequency uJc (right ordinate, thin dashed lines) as a function 
of the Prandtl number P in the case r] = 0.4 for the Coriolis 
numbers r = 5 x 10^, 10*, 1.5 x 10* and 10^ (from bottom 
to top). The thick lines correspond to expressions (15a) and 

(15c). 

In figure 2 expressions (15a,c) are compared with ac- 
curate numerical values which indicate that the general 



trend is well represented by expressions (15a, c). The 
same property holds for mg. For a rigorous asymptotic 
analysis in the case P = 1 including the radial depen- 
dence we refer to Jones et al. (2000). 




FIG. 3: Banana cells in a thin rotating spherical Buid shell 
cooled from within. Convection driven by centrifugal 
buoyancy is made visible by a suspension of tiny flakes 
which become aligned with the shear. (After Busse and 
Carrigan, 1976) 

It is evident from figure 2 that the agreement between 

expressions (15) and the numerical values deteriorates as 
low values of P are approached. This behavior is caused 
by the fact that instead of the thermal Rossby wave mode 
the inertial mode of convection becomes preferred at on- 
set for sufficiently low Prandtl numbers. It is character- 
ized by convection cells attached to the equatorial region 
of the outer boundary not unlike the "banana cells" seen 
in the narrow gap experiment of figure 3 . The equato- 
rially attached convection does indeed represent an in- 
ertial wave modified by the effects of viscous dissipation 
and thermal buoyancy. An analytical description of this 
type of convection can thus be attained through the in- 
troduction of viscous friction and buoyancy as pertur- 
bations as has been done by Zhang(1994) and by Busse 
and Simitev (2004) for stress-free and by Zhang (1995) 
for no-slip boundaries. According to Ardes et al. (1997) 
equatorially attached convection is preferred at onset for 
T < Ti where r; increases in proportion to P~^/^. 

A third form of convection is realized in the polar re- 
gions of the shell which comprise the two fiuid domains 
inside the tangent cylinder. Since gravity and rotation 
vectors arc nearly parallel in these regions (unless val- 
ues of 77 close to unity are used) convection resembles 
the kind realized in a horizontal layer heated from below 
and rotating about a vertical axis. Because the Coriolis 
force can not be balanced by the pressure gradient in this 
case, the onset of convection is delayed to higher values 
of R where convection outside the tangent cylinder has 
reached already high amplitudes. In the case of 77 = 0.4 
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FIG. 4: Convection in rotating spherical Euid shells in the 
cases T = 10*, J? = 4 x 10^, P = 1 (left column) and 

r = 5 X 10^, 7? = 8 X 10^ P = 20 (right column). Lines of 
constant Ur in the middle spherical surface, r = + 0.5 arc 
shown in the upper row. The plots of the middle row show 
streamlines, rd^/d(p = const., in the equatorial plane. The 

lowermost plots indicate lines of constant mean azimuthal 
velocity u^p in the left halves and isotherms of 6 in the right 

halves. 



the onset of convection in the polar regions typically oc- 
curs at Raylcigh numbers which exceed the critical values 
Rc for onset of convection outside the tangent cylinder 
by a factor of the order 4. Except for the case of very 
low Prandtl numbers the retrograde differential rotation 
in the polar regions generated by convection outside the 
tangent cylinder tends to facilitate polar convection by 
reducing the rotational constraint. A tendency towards 
an alignment of convection rolls with the North-South di- 
rection (Busse and Cuong 1977) can be noticed, but this 
property is superseded by instabilities of the Kiippers- 
Lortz type (for an experimental demonstration see Busse 
and Heikes (1980)) and by interactions with turbulent 
convection outside the tangent cylinder. 

Typical features of low and high Prandtl number con- 
vection are illustrated in figure 4. The columnar nature 
of convection does not vary much with P as is evident 
from the top two plots of the figure. At Prandtl numbers 
of the order unity or less, - but not in the case of inertial 




FIG. 5: The time- and azimuthally-averaged local Nusselt 
numbers at the inner spherical boundary Nui (left plot) and 

at the outer spherical boundary Nuo (right plot) as 
functions of the colatitude for P = 0.5, t = 1.5 x 10'' and 
R = 4,8, 10, 15, 20, 30, 25 x 10^ (from bottom to top at 90°). 



convection-, the convection columns tend to spiral away 
from the axis and thereby create a Reynolds stress which 
drives a strong geostrophic differential rotation as shown 
in the bottom left plot of the figure. This differential 
rotation in turn promotes the tilt and a feedback loop is 
thus created. At high values of P the Reynolds stress be- 
comes negligible and no significant tilt of the convection 
columns is apparent. In this case the differential rotation 
is generated in the form of a thermal wind caused by the 
latitudinal gradient of the axisymmetric component of 6. 

Among the properties of convection at finite ampli- 
tudes the heat transport is the most important one. Cus- 
tomarily its efficiency is measured by the Nusselt number 
which is defined as the heat transport in the presence of 
convection divided by the heat transport in the absence 
of motion. In the case of the spherical fluid shell two 
Nusselt numbers can be defined measuring the efficiency 
of convection at the inner and the outer boundary. 



Nui 



1 



P dB 

Ti dr 



Nuo = 1 • 



P dB 
r„ dr 



(16) 



where the double bar indicates the average over the 
spherical surface. In addition local Nusselt numbers 



Nui{6) = 1 



P_ dB 

n dr 



NuoiO) = 1 



p de 

rn dr 



(17) 

are of interest where only the azimuthal average is ap- 
plied, as indicated by the single bar, instead of the aver- 
age over the entire spherical surface. Examples of such 
measures of the dependence of the heat transport on lati- 
tude are shown in figure 5. This figure demonstrates that 
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at low supercritical Raylcigli numbers the heat transport 
occurs primarily across the equatorial region, but as R 
increases the heat transport in the polar regions takes 
off and soon exceeds that at low latitudes. This pro- 
cess is especially evident at the outer boundary. In the 
polar regions convection is better adjusted for carrying 
heat from the lower boundary to the upper one, and it is 
known from computations of the convective heat trans- 
port in horizontal layers rotating about a vertical axis 
that the value of Nu may exceed the value in a non- 
rotating layer at a given value of R in spite of the higher 
critical Rayleigh number in the former case (Somerville 
and Lipps, 1973). 

Outside the tangent cylinder the convective heat trans- 
port encounters unfavorable conditions in that the cylin- 
drical form of the convection eddies is not well adjusted 
to the spherical boundaries. This handicap is partly over- 
come through the onset of time dependence in the form 
of vacillations in which the convection columns expand 
and contract in the radial direction or vary in amplitude. 

At Prandtl numbers of the order unity and less an- 
other effect restricts the heat transport. The shear of the 
geostrophic differential rotation created by the Reynolds 
stresses of the convection columns severely inhibits the 
heat transport. To illustrate this effect we have plotted in 
figure 6 in addition to the Nusselt numbers the averages 
of the kinetic energy densities of the various components 
of the convection flow which are deflned by 

Ep = ^(1 V X (V$ X r) p), E, = 1(1 V§ x r f), 

(18a) 

Ep = ^(1 V X (V$ X r) Et = 1(1 V* X r 1^), 

(18b) 

where the angular brackets indicate the average over the 
fluid shell and $ refers to the azimuthally averaged com- 
ponent of while $ is defined by $ = $ — Analogous 
definitions hold for the energy densities of the magnetic 
field, 

Mp = ^(1 V X {Vh X r) Mt = 1(1 Vgxr \^), 

(19a) 

Mp = ^(1 V X {Vh X r) Mt = 1(1 Vgxr f). 

(19b) 

The total magnetic energy density M averaged over the 
fluid shell is thus given by M = Mp + Mi+ Mp + Mt. 

Figure 6 is instructive in that it demonstrates both, 
convection in the presence and in the absence of its dy- 
namo generated magnetic field. As is evident from the 
right part of figure 6 with M ~ relaxation oscillations 
occur in which convection sets in nearly periodically for 
short episodes once the energy Et of the differential ro- 
tation has decayed to a sufficiently low amplitude. But 
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FIG. 6: Decay of a dynamo for P = 0.5, t = 1.5 x 10*, 
i? = 1.2 X 10®, Pm = 0.5. The total magnetic energy density 

M (thick dotted Bnc, multipUod by a factor of 5), the kinetic 
energy densities Et (sohd Une), Ep (dotted line), Et (dashed 
line), and Nusselt rmmhers Nm (solid line) at the inner and 
Nuo (dashed line) at the outer spherical boundary (lower 
plot) are plotted as functions of time t. 



as soon as convection grows in amplitude, the differen- 
tial rotation grows as well and shears off the convection 
columns. After convection has stopped the differential 
rotation decays on the viscous time scale until the pro- 
cess repeats itself. 

As long as the magnetic field is present Et is suppressed 
owing to the action of the Lorentz force and high Nusselt 
numbers are obtained. The dynamo generated magnetic 
field thus acts in a fashion quite different from that of a 
homogeneous field which typically counteracts the effects 
of rotation and tends to minimize the critical value of the 
Rayleigh number when the Elsasser number. 



assumes the value 1 in the case of a plane layer (Chan- 
drasekhar 1961) or values of the same order in the case of 
a sphere (Fearn 1979) or in the related annulus problem 
(Busse 1983) . 
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FIG. 7: Convection driven dynamos as a function of the 
Rayleigh number R and the magnetic Prandtl number Pm 

for P = 0.1, r = 10^ (empty symbols) and for P = 1, 
T = 3 X 10'* (filled symbols). The symbols indicate chaotic 
dipolar (squares), hemispherical (triangles), and decaying 
dynamos (circles). 



5. CONVECTION DRIVEN DYNAMOS 

It appears that dynamos are generated by convection 
in rotating spherical shells for all parameter values as 
long as the magnetic Reynolds number is of the order 50 
or higher and the fluid is not too turbulent. Since Rm 
can be defined by Pm^/^E where the kinetic energy den- 
sity, E = Ep + Et+ Ep + Et, increases with R, increasing 
values of the Rayleigh number are required for dynamos 
as Pm decreases. In planetary cores the latter parameter 
may assume values as low as 10~^ or 10^^, but numeri- 
cal simulation have achieved so far only values somewhat 
below 10~^. The trend of increasing R with decreasing 
Pfn is evident in figure 7 where results are shown for two 
different values of t and P. It is also evident from the 
results in the upper left corner of the figure that an in- 
creasing R may be detrimental for dynamo action. This 
is a typical property of marginal dynamos at the bound- 
ary of dynamos in the parameter space (Christensen et al. 
1999, Simitev and Busse 2005). Since convection at the 
relevant values of R is chaotic this property also holds for 
the generated magnetic field. The magnetic energy must 
be finite near the boundary of dynamos in the parameter 
space because a sufficient amplitude of the fiuctuating 
components of the convection fiow can be obtained only 
when the magnetic field is strong enough to suppress the 
differential rotation. This property appears to hold even 
at high Prandtl numbers where differential rotation oc- 
curs only in the form of a relatively weak thermal wind. 

Even in its chaotic state convection at high values of 




p = 10, r = 5 X 10^, R = 5x 10^, = 10 




P = 5. r = 5 X 10^ R = 6 X 10^ P„ = 10 




p = 5, r = 5 X 103, P = 8 X 10^ P„ = 3 




FIG. 8: Various types of dynamo symmetry. The left 

column shows lines of Btp = const, (left half) and 
r sin Odeh = const, (right half). The right column shows 
surfaces of Br = const, at r = ro + 0.7. 



T exhibits a strong symmetry with respect to the equa- 
torial plane at high values of r, at least as long as con- 
vection in the polar regions is still weak. Because of this 
symmetry magnetic fields can be generated either with 
the same symmetry as the convection flow, i.e. h and 
g are antisymmetric with respect to equatorial plane in 
which case one speaks of a quadrupolar dynamo, or the 
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magnetic field exhibits the opposite symmetry with sym- 
metric functions h and g in which case one speaks of a 
dipolar dynamo. Of special interest are hemispherical 
dynamos (Grote and Bussc 2001) in which case the fields 
of dipolar and quadrupolar symmetry have nearly the 
same amplitude such that they cancel each other either 
in the Northern or the Southern hemisphere. Examples 
of typical dynamos with different symmetries are shown 
in figure 8. Here it is also evident that low P and high 
P dynamos differ in the structure of their mean toroidal 
fields. While their mean poloidal dipolar fields exhibit 
hardly any difference, the high P dynamo is character- 
ized by strong polar azimuthal flux tubes which are miss- 
ing in the low P case. The reason for this difference is 
that the radial gradient of the differential rotation in the 
polar regions is much stronger for high P than for low P. 

As the Rayleigh number increases and polar convec- 
tion becomes stronger the convection flow looses some of 
its equatorial symmetry and the magnetic field can no 
longer easily be classified. Usually the dipolar compo- 
nent becomes more dominant in cases of dynamos which 
started as either quadrupolar or hemispherical dynamos 
at lower values of R. 

It should be emphasized that quadrupolar and hemi- 
spherical dynamos arc usually oscillatory in that new az- 
imuthal flux tubes of alternating sign emerge in the equa- 
torial plane, move towards higher latitudes, while old flux 
tubes are dissipated in the polar region. There also exist 
dipolar dynamos which exhibit the same type of oscilla- 
tions. They are typically found in the parameter space 
near the region where hemispherical dynamos occur. To 
illustrate the dynamo oscillations additional plots sepa- 
rated at equal distances in time before the corresponding 
plot in the lower three rows of figure 8 are shown in fig- 
ure 9. The four plots in each case cover approximately 
half a period of oscillation. It must be realized, of course, 
that the oscillations are not strictly periodic since they 
occur in a turbulent system. In this respect they resem- 
ble the solar cycle with its 22 year period. Because the 
solar dynamo is believed to operate at the bottom of the 
solar convection the propagation of the dynamo wave is 
towards low latitudes on the sun. A remarkable feature 
of the dipole oscillation in figure 9 is that the polar flux 
tubes hardly change in time. The oscillation appears to 
be confined to the region outside the tangent cylinder. At 
a lower value of P„i the strong polar flux tubes even in- 
hibit the oscillation of the mean poloidal field as shown in 
figure 10. This case has been called the "invisibly" oscil- 
lating dynamo since at some distance from the boundary 
of the spherical fluid shell the oscillation of the dynamo 
can hardly be noticed. 

The relatively simple magnetic fields displayed in fig- 
ures 8, 9, and 10 should not be regarded as representative 
for high values of the Rayleigh number when R exceeds 
its critical value Rc by an order of magnitude or more. 
In that case the equatorial symmetry has nearly disap- 



peared owing to the dominant convection in the polar 
regions and non-axisymmetric components of the field 
tend to exceed axisymmetric ones as shown , for exam- 
ple, in figure 11. The strong growth with increasing Rm, 
of the non-axisymmetric magnetic flux is also evident in 
figure 12 when energies corresponding to different values 
of Pfn are compared. 

An important question is the average strength of the 
magnetic field generated by the dynamo process in de- 
pendence on the parameters. Two concepts are often 
used in attempts to answer this question. In astrophys- 
ical situations the equipartition between magnetic and 
kinetic energy is used as a guide. While such a balance 
may hold locally as, for example, in the neighborhood of 
sunspots, it is not likely to be applicable to global plan- 
etary magnetic fields. In the Earth's core, for example, 
the magnetic energy density exceeds the kinetic one by a 
factor of the order of 10^. The second concept is based 
on the property that an Elsasser number A of the order 
unity appears to correspond to optimal conditions for 
convection in rotating systems in the presence of an ap- 
plied nearly uniform magnetic field. While this idea may 
be useful as a first rough estimate, the results of numeri- 
cal dynamo simulations do not support this concept very 
well. The examples shown in figure 12 exhibit A- values 
differing by an order of magnitude and even larger vari- 
ations have been reported by Simitev and Busse (2005) 
and by Christensen and Aubert (2005). 

Another important feature demonstrated in figure 12 
is the change in the structure of the magnetic field with 
increasing Prandtl number. While for P = 5 and lower 
values the mean poloidal field is small in comparison with 
the fluctuating components, this situation changes dra- 
matically at about P = 8 for r = 5 • lO'"* such that for 
higher P the mean poloidal field becomes dominant. Usu- 
ally this field is dipolar. This change is associated with 
the transition from the geostrophic differential rotation 
generated by Reynolds stresses to the thermal wind type 
differential rotation caused by a latitudinal temperature 
gradient. While the magnetic energy M may exceed the 
total kinetic energy E by orders of magnitude in particu- 
lar for high Prandtl numbers, ohmic dissipation is usually 
found to be at most comparable to viscous dissipation in 
numerical simulations. But this may be due to the lim- 
ited numerically accessible parameter space. 

The brief introduction of this section to convection 
driven dynamos in rotating spherical fluid shells can only 
give an vague impression of the potential of numeri- 
cal simulations for the understanding of planetary mag- 
netism. Many more examples of such simulations can 
be found in the literature. Usually they have been mo- 
tivated by applications to the geodynamo and P = 1 
is assumed in most cases for simplicity. For some re- 
cent systematic investigations we refer to Christensen et 
al. (1999), Grote et al. (2000), Jones (2000), Grote 
and Busse (2001), Kutzner and Christensen (2000, 2002), 
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FIG. 9: Plots for the last three cases of figure 8 at the times n ■ At, n = 3, 2, 1, (left to right) before the times of figure 8 

with At =0.035,0.04,0.0025 (from top to bottom). 



Busse et al. (2003), Simitev and Busse (2005) and other 
papers referred to therein. 



6. APPLICATIONS TO PLANETARY DYNAMOS 
6.1 General Considerations 

In many respects it is too early to model the dynamo 

process in particular planets. The numerical simulations 
of the kind discussed in the preceding sections are still 
rather removed from the parameter regime relevant to 
planetary interiors. Only most recently attempts have 
been made to extrapolate results to high Rayleigh num- 
bers and high values of r (Christensen and Aubert, 2006). 
Moreover, only the most important physical parameters 
have been taken into account and relevant properties such 
as compressibility and other deviations from the Boussi- 
nesq approximation have been considered only in spe- 
cial models applied to the Earth's core (Glatzmaier and 



Roberts 1996) or to the sun (Brun et al. 2004). On the 
other hand, many essential parameter values of planetary 
cores are not sufficiently well known to provide a basis for 
the development of specific dynamo models. Much of the 
future progress of the field will thus depend on the mu- 
tual constraints derived from observational evidence and 
from theoretical conclusions in order to arrive at a better 
understanding of the workings of planetary dynamos. 

In table 2 some properties related to planetary mag- 
netism have been listed. The dipole moments of the plan- 
ets have been given as multiples of G ■ Rp where Rp de- 
notes the mean radius of the planet or satellite. Thus 
the numerical value indicates the field strength in Gauss 
in the equatorial plane of the dipole at the distance Rp 
from its center. Wc have included Mars in the table, 
although it does have an active dynamo at the present 
time. The strong magnetization of the Martian crust in- 
dicates, however, that a strong field dynamo must have 
operated in the early history of the planet. 

The most important question of the theory of plane- 
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FIG. 10: An "invisibly" oscillating dynamo with At = 0.04 in the case of P = = 5, r = 5 x 10^ and i? = 6 x The 
plots to the left show lines of constant in their left halves and meridional field lines, r sin ddgh = const., in their right 
halves. The plots to the right exhibits lines of constant Br on the surface r = ro + 0.4. The plots follow in clockwise sense 
starting with the upper left one such that they complete approximately a full cycle. 



Planet 
(Satellite) 


Planetary Radius 
Rp (10^ m) 


Dipole Moment 
(10-*TxJ^^) 


Core Radius 
(10® m) 


Angular Rate of 
Rot. n (IQ-^s^) 


Magn. Diffus. 
(mVs) 


Mercury 


2.439 


0.0025 


1.9 


0.124 


~ 2 


Earth 


6.371 


0.31 


3.48 


7.29 


2 


Mars 


3.389 


3+ 


~ 1.5 


7.08 


~ 2 


Jupiter 


69.95 


4.3 


56* 


17.6 


30* 


Saturn 


58.30 


0.21 


29* 


16.2 


4* 


Uranus 


25.36 


0.23 


18* 


10.1 


100* 


Neptune 


24.62 


0.14 


20* 


10.8 


100* 


Ganymede 


2.63 


0.0072 


0.7 


1.02 


4 



* Since the Giant Planets do not possess a well-defined boundary of a highly conducting core, 
values of the most likely region of dynamo activity are given. 
^ Remnant magnetism of the Martian crust appears to require an ancient dynamo 
with a field strength of at least 10 times the Earth's magnetic field. 



TABLE II: Planetary parameter values (after Stevenson (2003), Nellis et al. (1988), Connemey (1993) and various 

websites). 



tary magnetism is the dependence of the observed field 
strength on the properties of the planet. We have already 
discussed in the historical introduction various proposals 
for such dependences based on ad hoc assumptions, but 
these have been abandoned by and large. Only the con- 
cept of an Elsasser number A of the order unity is still 
frequently used. Stevenson (2003) shows that a value 
A = 0.3 fits most of the planets with a global magnetic 
field quite well, but in cases such as Mercury and Saturn 
only a value of the order 10~^ can be estimated for A. In 
order to save the concept different dynamo regimes must 
be assumed. Equilibration balances for weaker magnetic 



fields have been proposed as, for instance, in the case 
when the strength of convection is not suHicient to at- 
tain the A 1 balance. For details see Stevenson (1984). 

An important balance often invoked in discussions of 

planetary dynamos (see, for example. Stevenson (1979) 
and Jones (2003)) is the MAC-balance (Braginsky 1967) 
where it is assumed that Lorentz force, buoyancy force 
and pressure gradient arc all of the same order as the 
Coriolis force, while viscous friction and the momentum 
advection term are regarded as negligible in the equa- 
tions of motion. The ncglcction of the latter term is 
well justified for high Prandtl numbers, but it is doubt- 
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10* 



FIG. 11: Lines of constant Br at r = Vo + 1 for the 
dynamo with R = 2x 10®, r = 5 x 10^, P = = 1. 



ful for values of P of the order unity or less. Inspite 

of their smallncss in comparison with the Coriohs force, 
the divergence of the Reynolds stress can generate the 
most easily excitable mode of a rotating fluid, namely the 
geostrophic differential rotation which can not be driven 
by the Coriolis force, the pressure gradient or buoyancy. 
Of course, an excitation by the Lorentz force is possible 
in principle. The latter force, however, usually inhibits 
the geostrophic differential rotation. Since the mean az- 
imuthal component of the magnetic field is typically cre- 
ated through the shear of the differential rotation, the 
Lorentz force opposes the latter according to Lenz' rule. 

Christensen and Aubert (2006) have recently intro- 
duced a concept in which the equilibrium strength of 
planetary magnetism does no longer depend on the rate 
of rotation fl, but instead depends only on the avail- 
able power for driving the dynamo. They find that scal- 
ing laws can be obtained once Rayleigh number, Nusselt 
number and dimensionless buoyancy flux Q have been 
defined through quantities that no longer involve molec- 
ular diffusivities. This is a surprising result since it re- 
quires the presence of rotation and can not be achieved 
in a non-rotating system. The final estimate obtained for 
the strength of the magnetic field in the dynamo region 
fits the cases of the Earth and Jupiter quite well, but 
yields the result that Mercury's field could not be gen- 
erated by a buoyancy driven dynamo since the magnetic 
Reynolds number would be too small. In the cases of 
the outer planets relevant parameters are not sufficiently 
well known to draw definitive conclusions. 

Finally we like to draw attention to what has been 
called Stevenson's paradox, namely that a too high elec- 
trical conductivity is detrimental for dynamo action in 
planetary cores. According to the Wiedemann- Franz law 
of condensed matter physics the thermal conductivity k 
of a metal is proportional to its electrical conductivity 
times the temperature, k = aTL, where L has approxi- 
mately the same value for all metals. This is due to the 
fact that in metals electrons dominate both, the trans- 
ports of heat and of charge. As a consequence there is 
an upper bound on the electrical conductivity for which 
a dynamo driven by thermal convection is possible be- 
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FIG. 12: Kinetic Ex and magnetic Mx energy densities 
(upper part ) and viscous Vx and Ohmic Ox dissipations 

(lower part) as functions of Pm for convection driven 
dynamos for t = 5 x 10^, R= 10^ (left three columns), 
T = 3 X 10*, R = 3.5x 10® fright column) and Prandtl 
number as indicated in the boxes. The highest values of the 
Elsasser number for the cases P — 1, 5, 15, 2 arc A — 3.02, 
0.37, 0.49 and 1.12 respectively. The components Xp, Xt, 
Xp, Xt (where X = E, M, V, O) are represented by circles, 
squares, plus-signs and crosses, respectively. Kinetic energy 

densities and viscous dissipations are shown with light 
symbols, magnetic energy densities and Ohmic dissipations 

are shown with heavy symbols. The magnetic Prandtl 
numbers Pm = correspond to non-magnetic convection. 



cause a high value of a implies a low temperature change 

with increasing depth needed to carry the heat from the 
interior of the planetary core. But as soon as the temper- 
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FIG. 13: Sketch of Stevenson's paradox (after Stevenson 
2003). 



ature change with pressure falls below its isentropic value 
thermal convection disappears and a stably stratified en- 
vironment is obtained which makes it difficult to drive 
a dynamo by other types of motions such as composi- 
tional convection. Ignoring the latter possibility Steven- 
son (2003) has sketched the diagram shown in figure 13 
where the marginal nature of dynamos in terrestrial as 
well as in the "icy" planets is indicated. 

In the following brief characterizations of the mag- 
netisms of various planets and moons are given. In ter- 
restrial planets and satellites dynamos arc possible in the 
liquid parts of their iron cores. Such cores always include 
light elements which depress the temperature of freezing 
and can give rise to compositional buoyancy in the pres- 
ence of a solidifying inner part of the core. In Jupiter 
and especially in Saturn convection may also be driven in 
part by compositional buoyancy since hydrogen and he- 
lium are immiscible in certain regions and helium drops 
may rain out. In numerical dynamo models the different 
sources of buoyancy are usually not distinguished, how- 
ever. 



6.2 Mercury 

It has already become apparent that the origin of Mer- 
cury's magnetism is especially enigmatic. Besides the 
hypothesis of an active dynamo, thermoelectric currents 
and remnant magnetism of the crust compete as expla- 
nations. The latter possibility had been discarded for a 
long time since Runcorn's theorem states that a homoge- 
neous spherical shell can not be magnetized by an interior 
magnetic field in such a way that a dipolar field can be 



observed from the outside after the interior source has 
been removed (Runcorn, 1975). Aharonson et al. (2004) 
have shown, however, that plausible inhomogeneities of 
the crust could possibly explain a remnant magnetic field 
created by an ancient dynamo. 

The idea that a convection driven dynamo operates 
in a fluid outer core of Mercury has received some sup- 
port from observations of Mercury's librations (Margot 
et al., 2004) and from numerical simulations (Stanley et 
al., 2005). The amplitudes of the small periodic vari- 
ations (librations) of Mercury's rotation caused by the 
gravitational pull of the sun on the non-axisymmetric 
mass distribution of Mercury appear to be too large for 
a fully solidified planet. This suggests that the solid in- 
ner core is separated from the mantle by a fluid outer 
core in which convection flows may occur. The numeri- 
cal simulations of dynamos in the thin shell of the outer 
core have demonstrated that the dipole strength mea- 
sured from the outside may not be representative for the 
strength of the magnetic field at the core-mantle bound- 
ary and the criteria mentioned above may be satisfied 
after all. The question, however, whether there is a suf- 
ficiently strong source of buoyancy to sustain convection 
and the dynamo over the age of the planet has not yet 
been answered satisfactorily. 



6.3 Venus 

Space probes have clearly shown that there is no dy- 
namo operating presently in Venus. This result is sur- 
prising since Venus is very similar to the Earth in many 
other aspects. That the period of rotation is only 1/243 
of a day should not matter much since the Coriolis pa- 
rameter T is still huge. More important is the fact that 
the cooling of the planet is less efficient because of the 
apparent absence of mantle convection and plate tecton- 
ics. A solid inner core may not have yet started to grow 
in Venus and vigorous convection does not occur in the 
core, at least not at the present time. 



6.4 Mars and Moon 

Neither Moon nor Mars have an active dynamo, but 

the magnetized rocks in their crusts suggest that in their 
early history about 4 billion years ago dynamos may have 
been operating in the iron cores of these bodies. The 
small lunar core with a radius of the order of 350 km 
and the age of the magnetized lunar rocks put severe 
constraints on a possible dynamo origin of lunar mag- 
netism as discussed by Stcgman et al. (2003). In their 
model these authors try to accommodate in particular 
the apparent sudden onset of a lunar dynamo at a time 
3.5 Ga ago. Because of its larger iron core the possibil- 
ity of a temporarily operating Martian dynamo is more 
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likely, but the strong magnetization of parts of the Mar- 
tian crust requires a field of at least ten times the strength 
of the present geomagnetic field. For a recent review of 
the implications of Martian magnetism for the evolution 
of the planet we refer to Stevenson (2001). 

6.5 Jupiter 

The Jovian magnetic field is rather similar to the geo- 
magnetic one sharing with it a dipole axis that is inclined 
by an angle of the order of 10 degrees with respect to the 
axis of rotation. That the higher harmonics of the field 
are relatively stronger than those in the case of the Earth 
suggest that the Jovian dynamo is driven at a more shal- 
low layer than the Earth's outer core. This assumption 
is in accordance with the transition to metallic hydrogen 
which is expected to occur at a pressure of 140GPa cor- 
responding to a region with a radius of 0.9 of the Jovian 
radius (Nellis et al., 1996). The transition is not a sharp 
one as suggested by earlier models of Jupiter's interior, 
but a gradual one leading to a decrease of the magnetic 
diffusivity down to about 4m^/s. Since this value yields 
a rather high value for the Elsasser number A Steven- 
son (2003) has argued that the dynamo is most active 
above the region of highest conductivity. Jones (2003), 
on the other hand, accepts the value 20 for A, but must 
admit a high magnetic Reynolds number of the order 10^ 
at which an effective dynamo may no longer be possible. 
A detailed model of a convection driven Jovian dynamo 
should include effects of compressibility and depth de- 
pendent electrical conductivity, but none has yet been 
published. 



solely non-axisymmetric components are found in such 
situations. 

6.7 Uranus and Neptune 

A metallic liquid probably does not exist in Uranus and 
Neptune whose interiors consist mostly of ice and rocks. 
Because of the pressure dissociation of water deeper re- 
gions of these planets are characterized by an ionic con- 
ductivity which is lower than typical metallic conductiv- 
ities by between one and two orders of magnitude. The 
mixture of water, methane, ammonia and other ices is 
sufficiently fluid that convection can occur and that a 
dynamo is thus possible. The fact that the observed 
magnetic field do not show an alignment with the axis 
planetary rotation and that quadrupolar and octupolar 
components are comparable to the dipolar components 
of the fields has been attributed to a dynamo operating 
in a thin shell (Ruzmaikin and Starchenko, 1991; Stan- 
ley and Bloxham, 2004). This seems be hardly necessary, 
however, since even convection driven dynamos in thick 
shells often exhibit such fields when small scale compo- 
nents dominate as in the case of high Rayleigh numbers 
and low Prandtl numbers. See, for example, figure 11 
which shows a magnetic field rather similar to that of 
Uranus as displayed in the paper of Connerney (1993). 
A special difficulty for a convection driven dynamo in 
the interior of Uranus is caused by its low emission of 
heat. Holme and Bloxham (1996) suggest that a typical 
dynamo would involve more Ohmic dissipation than cor- 
responds to the heat flux emitted from the interior of the 
planet. 



6.6 Saturn 

Saturn's magnetic field with a strength of 0.2mT near 
the poles is much weaker than that of Jupiter which must 
be attributed primarily to its deeper origin in the planet. 
The transition pressure of 140GPa is reached only at 
about half the planetary radius. A property that has 
received miich attention is that the Saturnian magnetic 
field is nearly axisymmetric with respect to the axis of 
rotation. This property can not be interpreted as a con- 
tradiction to Cowling's (1934) theorem since it has long 
been shown that fields with arbitrarily small deviations 
from axisymmetry can be generated by the dynamo pro- 
cess (Braginsky 1976). Stevenson (1982) gave a reason- 
able explanation for the almost axisymmetric Saturnian 
field by demonstrating that the differential rotation in 
the stably stratified, but still electrically conducting gas 
shells above the dynamo region will tend to shear off all 
non-ax;isymmetric components of the field. The prob- 
lem may be a bit more intricate as has been pointed out 
by Love (2000) who showed that sometimes fields with 



6.8 Ganymede and other Satellites 

It came as a great surprise when the measurements 
of the Galileo spacecraft indicated that Jupiter's moon 
Ganymede possesses a global magnetic field for which an 
active dynamo inside the satellite seems to be the only 
realistic explanation (Kivelson et al., 1996), although the 
possibility of a remnant magnetism can not easily be ex- 
cluded (Crary and Bagenal 1998). It must be kept in 
mind that Ganymede is the largest satellite in the solar 
system which exceeds even Mercury in size. Neverthe- 
less the estimated radius of its iron core is only about 
660 km and it is hard believe that it could still be partly 
molten unless it contain a lot of radioactive potassium 40 
or Ganymede was captured into a resonance in its more 
recent history (Showman et al. 1997). The fact that 
Ganymede's magnetic moment is nearly aligned with the 
ambient Jovian magnetic has led to the suggestion that 
the ambient field could aid significantly Ganymede's dy- 
namo. Sarson et al. (1997) have investigated this ques- 
tion with the result that the ambient field is too weak to 
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exert much influence. On the other hand, in the case of To 
which is much closer to Jupiter the interaction between 
the ambient field and the liquid iron core can explain 
the observed structure of the magnetic field without the 
assumption of an active dynamo. 

There are no other satellites in the solar system where 
an active dynamo must be suspected. The siibstantial 
magnetization of many iron meteorites suggest, however, 
that several differentiated proto-planets have had dy- 
namos in the early days of the solar system. 

7. CONCLUDING REMARKS 

It is apparent from the above discussions that dynamo 
theory does not yet have much specific information to 
contribute to the interpretation of the observed magnetic 
properties of planets and satellites in the solar system. 
Even possible interactions between thermal and compo- 
sitional buoyancies have not yet been taken fully into 
account (Glatzmaier and Roberts 1996, Busse 2002b). 
Some typical properties are already apparent, however, 
as for instance: 

• Dynamos that are dominated by a nearly axial 
dipole and exhibit a magnetic energy that exceeds 
the kinetic one by orders of magnitude as in the case 
of the Earth are typical for high effective Prandtl 
numbers as must be expected when convection is 
primarily driven by compositional buoyancy. The 
rather Earth like appearance of the dynamo of 
Glatzmaier and Roberts (1995) is in part due to 
the high value of P used in their computational 
model. 

• Dynamos exhibiting strong higher harmonics are 
more likely driven by thermal convection corre- 
sponding to Prandtl numbers of the order unity or 
less. 

• Dynamo oscillations are a likely phenomenon in 
the presence of a sufficiently strong differential ro- 
tation. They may not always be apparent in the 
poloidal magnetic field seen at a distance from the 
dynamo region. 

• Considerations based on the effects on convection 
of an imposed homogeneous magnetic field can not 
directly be applied to the case of convection driven 
dynamos. 

More details on the parameter dependence of dynamos 
will certainly emerge in the future as the increasingly 
available computer capacity will allow extensions of 
the parameter space accessible to computer simulations. 
Space probes such as MESSENGER in the case of Mer- 
cury will provide much needed detailed information on 
planetary magnetic fields. Eventually we may learn more 



about their variation in time which is one of the most in- 
teresting properties of planetary magnetic fields. 

As in the case of stellar magnetism where the study of 
star spots and stellar magnetic cycles is contributing to 
the understanding of solar magnetism it may eventually 
become possible to learn about the magnetism of extra- 
solar planets and apply this knowledge for an improved 
understanding of solar system dynamos. Undoubtedly 
the field of planetary magnetism will continue to be an 
exciting one! 
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